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probe D5-branes in global Anti-de Sitter space. The dual field theory is A/" = 4 theory 
on IR X S''^ with added fiavour. The branes undergo a geometrical phase transition in 
this geometry as function of the bare quark mass in units of with R the 
radius. The meson spectra are obtained from fiuctuations of the brane probes. First, we 
study them numerically for finite quark mass through the phase transition. Moreover, 
at zero quark mass we calculate the meson spectra analytically both in supergravity 
and in free field theory on IR x S''^ and find that the results match: For the chiral 
primaries, the lowest level is given by the zero point energy or by the scaling dimension 
of the operator corresponding to the fiuctuations, respectively. The higher levels are 
equidistant. Similar results apply to the descendents. Our results confirm the physical 
interpretation that the mesons cannot pair-produce any further when their zero-point 
energy exceeds their binding energy. 
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1. Introduction 



The AdS/CFT correspondence in its original form involves the near-horizon limit of 
D3-branes. This limit leads to the AdS^ x geometry, where the Anti-de Sitter factor 
of the geometry is realized as the Poincare patch of AdS^. 

Following [1], there have been extensive studies of probe D7-branes in this geometry 
which on the dual field theory side leads to added hypermultiplets in the fundamental 
representation of the gauge group. For preserving J\f = 2 super symmetry, the probe 
D7-brane has to wrap a subspace which is asymptotically AdS^ x 5*^ near the boundary. 
The complete meson spectra arising from the fluctuations of the probe brane, ie. for 
all scalar, fermion and vector modes, have been found in [2]. For instance, for scalar 
mesons the result is [2] 

or T m 

M,{nJ) = —^{n + l + l){n + l + 2), — = V27r-^, (1.1) 

where L is the embedding coordinate of the D7 brane, R is the AdS radius, rUg is 
the quark mass and A the 't Hooft coupling, n is the main quantum number and / is 
the quantum number of the SO {A) symmetry associated with the asymptotically 
wrapped by the D7-brane probe near the boundary. The spectrum has a large degener- 
acy since it depends only on the combination (n + l). This is expected from the M = 2 
supersymmetry of the system. For general Dp/Dq-brane systems, similar spectra were 
found in ref. [3, 4]. 

By embedding a D7-brane probe into deformed versions of AdS^ x 5*^, physical 
phenomena may be described. Examples are chiral symmetry breaking, which is ob- 
tained, by embedding a D7-probe into a confining background [5], and a first-order 
meson-melting phase transition at finite temperature which arises when embedding 
a probe brane into the AdS-Schwarzschild black hole geometry [5, 6, 7, 8]. For the 
AdS-Schwarzschild geometry, a topology-changing phase transition occurs, depending 
on the ratio of quark mass over temperature, between branes that reach the black hole 
horizon and those that do not. 

Supersymmetric embeddings of D5-brane probes wrapping an AdS4^ x S'^ of AdS^ x 
have first been studied in [9]. In this case, the dual field theory is a superconformal 
'defect' theory [10] in which the additional hypermultiplets are confined to a (2+1)- 
dimensional subspace. 

In [11] (see also [12]), the authors investigate the embedding of D7- and D5-brane 
probes into global AdS. Global AdS is dual to A/" = 4 Super Yang-Mills theory on 
R X S*^. Moreover these authors study the thermodynamical properties of these systems 
by considering brane probes in global thermal AdS, where also the time direction is 
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compactified, such that the field theory is defined on x S^. Again in these geometries 
they find a topology-changing phase transition, depending on the ratio of quark mass 
over radius, between those embeddings that reach the and those that do not. 
Using critical exponents, the authors show (see also refs. [13, 14, 15]) that the phase 
transition is third order for D7-brane probes, while it is first order for D5-brane probes. 
They also consider fiuctuations of the brane and calculate the mass of the lowest-lying 
scalar meson mode. The authors find that the spectrum has a kink at the phase 
transition. Moreover, by making use of the Polyakov loop, they study the physical 
properties of the low-energy phase. In this phase, the zero-point energy of the mesons 
- which is due to the finite volume of the 5"^ - is larger than their binding energy. 
This means in particular that the mesons are deconfined in the sense that they cannot 
pair-produce any more. String breaking is no longer possible. 

In this paper we investigate the mesonic spectrum on R x S*^ in further detail. We 
look at scalar and vector fiuctuations for both the D7- and the D5-brane probe cases. 
For the D7 brane probe, we establish the complete bosonic spectrum on the gravity 
side. Our most important result concerns the spectrum at vanishing quark mass, for 
which we perform analytical calculations both on the gravity and on the field theory 
side. On the gravity side, we map the fiuctuation equations of motion to equations of 
Schrodinger type. For the chiral primaries, we find that the energy of the ground state is 
given by the dimension of the operator dual to the fiuctuations. The higher fiuctuation 
modes, labelled by the quantum number n, are equidistant. For the descendents, we 
also find an equidistant spectrum. However the energy of the ground state is no longer 
equal to the dimension of the operator. 

We then turn to the field theory side, and in particular to the free M = 2 theory 
in 3 + 1 dimensions obtained by adding a fiavour hypermultiplet to the original A/" = 4 
theory. We expand the fields on the within R x S*^. The spectrum of the mesonic 
composite operators is obtained by combining the expansions of the component fields, 
making use of the appropriate Clebsch-Gordan coefficients. This requires a careful 
analysis of the transformation properties of the mesonic composite operators under the 
antipodal map.^ For both the chiral primaries and the descendents, we find the same 
result for the spectrum as in the gravity calculation. As an example, let us quote our 
result for scalar D7-brane mesons, which is 

M{n,lJ) = h?> + 2n + l + l). (1.2) 
K 

Here, R is the radius of the in the field theory directions, which gives rise to the 

^The antipodal map A : 5*" — > S"", defined by A{x) = —a;, sends every point on the sphere to its 
antipodal point. 
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conformal mass 1/ R. Note that since this mass arises from the background metric rather 
than from the D7-brane boundary conditions, (1.2) is independent of the 't Hooft 
couphng A, while (1-1) is (9(1/a/A). Moreover, (1.2) depends on the main quantum 
number n, as well as on the quantum numbers / and /, where I refers to the internal 
S*^, asymptotically wrapped by the D7-brane, while / refers to the S*^ in IR x in 
the field theory directions. The lowest mode, ie. the zero-point energy, is given by the 
conformal dimension of the dual operator, A = 3. 

It is remarkable that the free field calculation and the gravity calculation of the 
meson spectrum agree. This is of course due to non-renormalization theorems which 
hold also in = 2 theory. From the physics perspective this confirms the physical 
interpretation, already advocated in [11], that the mesons cannot pair-produce even at 
strong coupling when confined into a small volume such that their zero-point energy is 
larger than their binding energy. 

We note that when setting / to zero, our new result (1.2), which in this case 
depends on the combination 2n + /, is less degenerate than (1.1) which depends on 
n + I. This linked to the fact that for supersymmetric field theories on S*^, scalars, 
fermions and vectors in the same multiplet have different conformal mass. Also, the 
supersymmetry transformations of the field theory fermions are modified to contain 
curvature-dependent terms. For A/" = 4 Super- Yang-Mills theory on IR x S^, these 
issues have been studied in detail in [16]. We expect similar results to hold for the 
M = 2 theory considered here. 

Our calculation shows - in a simple example - that it is possible to directly compare 
gravity and field theory calculations in top-down approaches to holographic models of 
physical relevance. Such top-down models have been used widely recently to holograph- 
ically describe physical phenomena such as superfluidity, quantum critical points and 
transport properties [17, 18, 19, 20]. These models have the advantage - as compared 
to bottom-up models - that the dual field theory is explicitly known. We hope that 
further field-theory studies will follow soon for a more detailed comparison between the 
weak and strong coupling aspects of a given model. 

This paper is organized as follows. In section 2 we introduce the general setup of 
global AdS and its probe brane embeddings, and present the phase transitions for the 
D7 and D5 brane cases. In section 3 we study the probe brane fluctuations as function 
of the bare quark mass and discuss the behaviour of the meson spectrum for both the 
phase transition and the case of vanishing quark mass. In section 4, we analytically 
compute the meson spectra at zero quark mass on the gravity side. We obtain the full 
bosonic spectrum for the D7 brane fluctuations, as well as some characteristic examples 
for the D5 brane case. In section 5 we present the field theory calculation at zero quark 
mass and show that it agrees with the gravity results. We end with concluding remarks 
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in section 6. 

2. General Setup 

Our starting point is AdSs x in global coordinates: 

ds^ = -(1 + w''/R^)de + W^dnl + ^ ^^2/^2 + ^'^^5 • (2-1) 

It is convenient to introduce the following radial coordinate: 

u = ]^{w + ^WTw^) . (2.2) 
The metric (2.1) in these coordinates is given by: 




Note that u > R/2. In this way the transverse has a ball of radius R/2 sitting at 
the origin. 

2.1 Probe D7— brane. 

Let us introduce a probe D7-brane. To this end it is convenient to write the metric of 
unit in the following coordinates: 

ds% = de^ + cos^ ed^l + sin^ Odcf^ . (2.4) 

Now if we let the D7-brane be extended along the AdS^ part of the geometry and wrap 
an S'^ <Z , the radial part of the corresponding DBI lagrangian is given by: 

^ ^ " i£) - ' ^ v^TT^I^ . (2.5) 

Possible embeddings split into two classes Minkowski embeddings and "ball" embed- 
dings correspondingly wrapping shrinking 5*^ cycles in the and AdSs parts of the 
background [12] (look at figure 1). The two classes are separated by a critical em- 
bedding which has conical singularity at the ball (represented by the dashed curve in 
figure 1). According to the AdS/CFT dictionary the bulk dynamics of the scalar 9{u) 
encodes the dynamics of the dual gauge invariant operator. In particular one can read 
off the source and the vev of the operator from the asymptotic behaviour of 9{u). In 
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Figure 1: Blue curves correspond to Minkowski embeddings and Red curves correspond to 
"ball" embeddings. The black dashed line corresponds to the critical embedding. 



our case the operator is the fundamental bilinear and the source and vev of the operator 
correspond to the mass of the hypermultiplet and the fundamental condensate. The 
precise dictionary has been derived in ref. [21], where an elegant renormalization pre- 
scription has been offered. Let us briefly review the results of refs. [21, 12] in a slightly 
modified form relevant to our notations. For large u the solution to the equation of 
motion derived from (2.5) has the following expansion: 

e{u) = + —y- - -— log - + . . . . (2.6) 

After following the renormalization prescription outlined in ref. [21] the condensate of 
the theory can be calculated in terms of the parameters 0o,62- For our choice of radial 
coordinate the result is given by: 

(V^V') oc -2^2 + 1 + ^0 log ^o; • (2.7) 

If we split the space corresponding to du^ + u'^dfl'^ to x and define radial 
coordinates p = u cos 9 and L = u sin 9 the DBI lagrangian in this coordinates is given 
by: 
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The profile of tlie D7-brane embedding L{p) lias the following asymptotic behaviour 
at large p: 

Lip)=m+-^-—\ogp + ... . (2.9) 

Further more one has the relations m = 6qR and ci = R^{62 — ^0/6). The condensate 
of the theory is then given by: 

(ipij) oc -2ci + mi?^ log(m/i?) = -2c , (2.10) 

where we have defined a new parameter c proportional to the condensate. Note also 
that the bare mass of the hypermultiplet is proportional to m the exact relation is 
rriq = m/2TTa' . 

After solving numerically for the D7-brane embeddings we can generate a plot of 
the equation of state c(m) presented in figure 2, where we have used dimensionless 
parameters fh = m/R and c = c/R^. As one can see there is no apparent multi- valued 
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Figure 2: A plot of the condensate — c versus the bare mass m. The states corresponding to 
"ball" embeddings are presented by the red curve and the states corresponding to Minkowski 
embeddings are presented by the blue curve. 

region near the transition from Minkowski to "ball" embeddings. In fact it has been 
shown that the phase transition is continuous and is of a third order [11]. 

Let us comment briefly on the physical meaning of the parameter m = m/ R and 
in particular its R dependence. Consider a constant u = A ^ R slice of the AdSs x 5*^ 
space-time. The induced metric has the asymptotic form: 

A^ 

dsl = —.{dt^ + R^dn^) + R^dnl . (2.11) 
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Here A is specifying the energy scale of the dual gauge theory. It is natural to identify 
the radius of the AdSs space-time R with the radius of the three sphere i?3 where the 
holographically dual field theory is defined. Furthermore we know that m = 27Ta'mq, 
where rriq is the bare mass of the hypermultiplet. However the quantity mq/R^ is 
obviously not a dimensionless parameter. To clarify this let us consider a rescaling 
A — 7- 7A of the energy scale. The metric in equation (2.11) becomes: 

A^ 

ds^A = -^h^dt^ + 'j^R^dnl) + R^dnl . (2.12) 
R 

We can always scale the time coordinate t — )■ t/'-f to compensate for the rescaling of 
the energy scale. Note that this suggests that the bare mass of the hypermultiplet 
should also be rescaled rriq — 'yrriq. Furthermore the radius of the three sphere is now 
i?3 = '-fR and hence we can write: 7 = R3/R. For the parameter m we obtain: 

_ 7m,(27raO _ 27ra' _ it niqR^ 

m - mqRs^ - -^-^ , (2.13) 

where we have used that R^ is related to the t'Hooft coupling of the dual gauge theory 
via: R^ = 2Xa'. Equation (2.13) specifies the physical meaning of the parameter rh. 
Now we can interpret the phase transition as taking place at constant bare mass 
and varying radius of S^. Note that varying R^ corresponds to varying the Casimir 
energy of the dual gauge theory and hence it is essentially a quantum phase transition. 



2.2 Probe D5— brane 

Let us review the case of a probe D5-brane studied in ref. [11]. To this end it is 
convenient to write the part of the geometry in the following coordinates: 

dsls = dip^ + cos^ ipdnl + sin^ tpdClj . (2.14) 

Next we define: r = ucosip and / = usinip. The metric (2.3) can then be written as: 

ds' = -^,{l + ^y de + ^ (1 - ^) ' [d-' + sin^ a(rf/3^ + sin^ /^^T^)] 

+ ^ [dr^ + r^dnl + di^ + fdn''^ . (2.15) 

Now if we let the D5-brane be extended along the t, /3, 7, r, Q2 directions and has a 
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non-trivial profile along /, the corresponding DBI lagrangian is given by: 



£ oc 1 + 



1 



rVl + / 



12 



(2.16) 



It is easy to show that the solution to the equation of motion has the following asymp- 
totic behaviour at large r: 



l{r) 



m H h • . . . 

r 



(2.17) 



Note that unlike the D7-brane case there is no extra logarithmic term [21] in the 
asymptotic expansion of /(r). One can then directly relate the coefficient c to the fun- 
damental condensate of the theory {{qq) cx — c). It is convenient to define dimensionless 
quantities c = c/ and m = m/R. 

After solving numerically the equation of motion of the probe D5-brane we obtain 
the plot of — c versus fh presented in figure 3. As one can see there is a multi-valued 
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Figure 3: A plot of the condensate — c versus the bare mass fh for the D5-brane probe. 
The states corresponding to "ball" embeddings are presented by red curves and the states 
corresponding to Minkowski embeddings are presented by blue curves. The zoomed in plot 
suggests a first order phase transition in the dual gauge theory. 



- 9 - 



region near the state corresponding to the critical embedding. This suggests that in 
this case the corresponding phase transition in the dual gauge theory is a first order 
one. The value of the critical bare mass at which the phase transition takes place 
can be obtained using the equal area law as illustrated in the second plot in figure 3. 
The authors of ref. [11] showed that the order of the phase transition can be inferred 
analytically, by calculating appropriate critical exponents. 

3. Meson spectra 

Some aspects of the spectrum corresponding to fluctuations of the D7-brane embed- 
dings were studied in ref. [11], where the lowest mode for each of the fluctuations was 
considered. In particular, the authors considered the scalar meson frequency squared 
corresponding to the fluctuation along the profile of the embedding (the coordinate L 
in our case) and found that as function of the quark mass parameter m, it has a kink 
precisely at the critical value where the phase transition occurs. This kink was 
shown to be consistent with the fact that the phase transition is third order. On the 
other hand, the spectrum of fluctuations along the coordinate (p was shown to be a 
smooth function of rh across the phase transition. Moreover, the authors of ref. [11] 
pointed out that for large rh ^ 1 (ie. for large 5*^ radius) the spectrum matches with 
the meson spectrum of the D3/D7 brane intersection obtained in ref. [2]. 

The authors of ref. [1 1] explain the behaviour of the lowest meson mode as follows: 
Consider the case that at fixed quark mass, the system is squeezed into smaller and 
smaller volume. When the compactification radius R passes through its critical value 
R*, the mesons are expected to deconfine as their finite volume zero point energy within 
becomes larger than their binding energy. Beyond the phase transition the theory is in 
a deconfined phase, in the sense that it cannot pair-produce. 

Below we provide a more detailed analysis of the spectrum of meson-like excitations 
by considering higher excited modes (n > in the notations of ref. [2]). Our study 
shows that at the phase transition, the qualitative behaviour of the higher modes is the 
same as that of the ground state, as far as the appearance of the kink is concerned. For 
values m < m*, we continue to find a discrete spectrum. In particular, for the limit of 
zero bare mass (m = 0), we are able to perform an analytic calculation, presented in 
detail in section 4 and 5, mapping the fluctuation equation of motion to an equation 
of Schrodinger type, and find that the discrete spectrum is equidistant. Moreover, for 
the fluctuations dual to chiral primaries we find that the energy of the ground state 
is given by the conformal dimension of the dual operator (in units of 1/R). A similar 
structure emerges for the fluctuations of a D5 brane probe. 



- 10 - 



Moreover, for the D5 brane probe fluctuations we also investigate the m dependence 
and show that the quahtative behaviour near the critical embedding is consistent with 
the phase transition being first order. In particular, the spectrum develops tachyonic 
states exactly at those points moo of the multi-valued region of the condensate c(m) 
where the slope of the condensate diverges. 

We also turn on the Kaluza-Klein modes of the D7-brane in the part of the 
geometry corresponding to operators with non-vanishing R-charge. We find that the 
large degeneracy of the spectrum of the D3/D7 system considered in ref. [2] is lifted for 
the theory on S^. In the limit of large rfi (large radius of S^) this degeneracy is restored. 
In the limit of zero bare mass the degeneracy is only partially restored. Furthermore, 
once again the energy of the ground state is determined by the engineering dimension 
of the operator (A = 3 + /) and the spectrum is equidistant. 

The outline of our calculations and results in this section is as follows: We first 
present our numerical results for fluctuations of a probe D7 — brane along the coordi- 
nate L. We describe the qualitative behaviour across the phase transition and discuss 
the structure of the spectrum at vanishing bare mass. Next we present our numerical 
results for the spectrum of fluctuations along for the D7-brane case, discussing also 
the Kaluza-Klein modes. In the limit of vanishing bare quark mass, we compare our 
numerical results to the analytical results obtained below in section 4 and 5. Moreover 
we analyze the spectrum of fluctuations of the D5-brane. We discuss both fluctua- 
tions along the / coordinate and along the Neumann-Dirichlet coordinate (the polar 
coordinate of C AdS^). 

3.1 Fluctuations of the D7— brane embedding. 
3.1.1 Fluctuations of the transverse scalars. 

In order to study the light meson spectrum, we look for the quadratic fluctuations of 
the D7-brane embedding along the transverse directions parametrized by L, 0. To this 
end we expand: 

L = L + 27ia'6L; = , (3.1) 

in the lagrangian (2.8) and leave only terms of order (27rQ;')^. After some calculations 
we obtain the following lagrangian for the quadratic fluctuations along L: 



2 1 1 



dp yi + L'"^ 



SL\ (3.2) 



where Gll is the corresponding component of the metric of the AdS^ x background 
and 

gap ig 

the inverse of the induced metric on the worldvolume of the D7-brane. 
Next we obtain the equation of motion for SL and consider an ansatz 6L = e*'^*/i(p). 
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Finally we solve the resulting differential equation for h{p) numerically and require 
normalizability of the solution, which leads to a quantized spectrum. The corresponding 
spectrum is presented in figure 4, where the dimensionless quantity Co = uR has been 
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Figure 4: A plot of tfie meson mass versus tfie bare quark mass m for the spectrum 
of D7 fluctuations along L. The black dashed lines correspond to the spectrum of D3/D7 
intersection [2] given by equation (3.3), one can observe the good agreement at large m. At 
the critical bare mass (represented by the vertical dashed line) the spectrum has a kink. 
Interestingly at m = the spectrum is discrete and equidistant taking only integer values 
given by equation (3.4). 

defined. At large values of fh (large radius of S^) the spectrum approaches the spectrum 
of the D3/D7 intersection given by 

^ = ^Vin + l)in + 2). (3.3) 

At the critical bare mass (represented by the vertical dashed line in figure 4) the 
spectrum has a kink as pointed out in ref. [11]. Below the phase transition the spectrum 
has a discrete quasi-equidistant structure which becomes exact at vanishing bare mass: 

u = {2n + 3)^. (3.4) 

Note that in (3.4), the energy of the ground state {n = 0) is given by the engineering 
dimension A = 3 of the operator dual to the fluctuations considered. 

As the authors of ref. [11] argue in this phase the mesons are deconfined. The 
equidistant structure of the spectrum then arises from the Kaluza-Klein tower of the 
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collection of fields on corresponding to the field content of the deconfined mesons. 
We refer the reader to sections 4 and 5 for analytic derivations of equation (3.4) both 
on the supergravity and on the field theory side of the correspondence. 

Note that the dependence of the spectrum on the 't Hooft coupling A is through 
the parameter m. According to equation (3.3), at large bare mass m the dimensionless 
quantity u depends linearly on the parameter m and hence in view of equation (2.13) 
is proportional to In agreement with figure 4, below the phase transition, i.e. 

at small values of m, the spectrum changes weakly with m. Asymptotically, for m — )■ 
the spectrum becomes independent of rh and hence of A. 

Let us now focus on the spectrum of fluctuations along (p. After expanding to second 
order in a' in the DBI lagrangian (2.8) and substituting the ansatz 5(j) = e''^%m{S^) f (p) 
into the resulting equation of motion we obtain: 




/(/ + 2) 

p2 



/(p) = . (3.5) 



Note that the quantum number / represents Kaluza-Klein modes on 5*^ C and labels 
representations of the global SU(2) R-symmetry of the dual gauge theory. Let us first 
focus on the spherically symmetric case (/ = 0). Solving equation (3.5) numerically and 
requiring normalizability of the solution, we obtain the meson mass u versus the bare 
quark mass m presented in figure 5. We see that for large bare mass (or equivalently for 
large radius of S^), the spectrum matches the spectrum of the D3/D7 system described 
by equation (3.3). In the limit of zero bare mass (m = 0), the spectrum is equidistant 
and is again identical to spectrum of fluctuations along L described by equation (3.4). 
At the phase transition (the vertical dashed line in figure 5) the spectrum is a smooth 
function of rh. 

Our next goal is to study the spectrum at non-zero /. Note that in the "flat 
case" considered in ref. [2], ie. with a Minkowski boundary in Poincare coordinates, the 
spectrum has degeneracy in n and / and is described by generalization of equation (3.3): 

2t/7 / 

io = ^Vin + l + l){n + l + 2) . (3.6) 

On the other hand, the theory on , which we consider here, is not supersymmetric 
for generic values of in and we do not expect the spectrum to be degenerate. Indeed 
our numerical investigations confirm this. As an illustration, in figure 6 we provide the 
spectrum of the ground state (n = 0) for various values of /. As one may expect, 
at large ffi the spectrum again matches the spectrum of the D3/D7 system given by 
equation (3.6). We also compared to the / = 0, n 7^ result in figure 5 and verified that 
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Figure 5: A plot of the meson mass w versus the bare quark mass rh for the spectrum of D7 
fluctuations along ip. At large fh the spectrum is described by equation (3.3). At zero bare 
mass (tti = 0) the spectrum is equidistant and is described by equation (3.4). Note that the 
spectrum is a smooth function of m across the phase transition. 
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Figure 6: D7 brane fluctuations: A plot of the meson mass w versus the bare quark mass 
m for various values of I. At large m the spectrum is given by equation (3.6). At zero bare 
mass the spectrun is described by a; = 3 + Z. 

for generic values of m the spectrum is not degenerate. This is particularly evident in 
the deconfined phase (the red colored curves). Interestingly, the spectrum at m = is 
described by cD = 3 + /, which suggests that the spectrum again has some degeneracy. 
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This fits with the fact that supersymmetry is restored at m = 0. In fact, by taking 
the limit L — )■ in equation (3.5) and analyzing the resulting equation of motion 
analytically, we find that the spectrum at m is given by: 

w = (3 + 2n + /)4 • (3.7) 

We refer the reader to section 4 for a detailed derivation of the spectrum at m = 0. 
Notice however that the combination 3 + / in equation (3.7) is equal to the engineering 
dimension of the corresponding dual field theory operator. Indeed if we look at the 
asymptotic form of the equation of motion (3.5) for large p we obtain: 

f"{p) + -f\p) - ^^^f{p) = . (3.8) 
P P^ 

The solution of equation (3.8) is of the general form: 

/(p) = p'<5m + p-('+2)5c (3.9) 

and hence according to the standard AdS/CFT dictionary, the corresponding gauge 
invariant operator indeed has conformal dimension A = 3 + /. This suggests that for 
generic meson-like excitations the spectrum at zero bare mass is given by: 

w = (A + 2n)4, (3.10) 

where A is the conformal dimension of the field operator dual to the lowest fluctuation 
mode in agreement with the standard field/operator dictionary. Note also the increase 
in units of 2n between the equidistant levels. We study the spectrum in more detail in 
sections 4 and 5 below. For completeness and to check the validity of equation (3.10), 
we now turn to analyzing the spectrum of fluctuations of a probe D5-brane in the next 
subsection. 

3.2 Fluctuations of the D5— brane embedding. 

In order to study the meson spectrum corresponding to fluctuations of the D5-brane 
along the coordinate /(r), defined above (2.15), we expand /(r) = l{r) + (27ra')(5/(r) 
in the DBI Lagrangian (2.16) and leave only terms of order a'^. Next we obtain the 
corresponding equation of motion and consider an ansatz 5l{r) = e^'^^ri{r). After solving 
the equation of motion numerically and imposing regularity of the solution, we obtain 
the meson spectrum as a function of the bare quark mass. Our results for the ground 
state and the first few excited states are presented in figure 7, where u = uR as before. 
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Similarly to the D7-brane case, the spectrum for large m matches the spectrum of the 



n=2 





Figure 7: A plot of the meson mass u} versus the bare quark mass rh for the spectrum of D5 
brane fluctuations corresponding to l{r). At large m the spectrum is described by equation 
(3.11). At zero bare mass {rh = 0) the spectrum is equidistant and is described by equation 
(3.10) with A = 2. Near the phase transition the spectrum becomes tachyonic at the points 
where the slope of the — c versus rh plot from figure 3 diverges. 

D3/D5 intersection studied in refs. [3, 4] where the following analytic expression has 
been obtained: 

u^ = ^Vin + l/2)in + 3/2). (3.11) 

For n = 0,1, 2, (3.11) is represented by the dashed black lines in figure 7. Interestingly 
at vanishing bare mass (m = 0), the spectrum is discrete and equidistant and is de- 
scribed by equation (3.10) with A = 2. This is to be expected since fluctuations along / 
correspond to fluctuations of the fundamental condensate of the dual gauge theory and 
for the defect-field theory that we consider (2+1 dimensional defect), the condensate 
operator is a dimension two operator (A = 2). We refer the reader to sections 4 and 5 
below for an analytic derivation of the spectrum at zero bare mass both on the gravity 
and on the field theory side of the correspondence. 
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The second plot in figure 7 represents the structure of the spectrum near the phase 
transition. From the first order phase transition pattern of the plot of — c versus rh 
presented in figure 3, one may expect that the theory becomes unstable at the points 
where the slope of the fundamental condensate diverges. As we can see from the graph 
for the ground state {n = 0) in figure 7, the spectrum drops to zero exactly at these 
points. We have verified numerically that beyond these points the spectrum of the 
ground state is indeed tachyonic. 

For completeness of our study and to verify the stability of the D5-brane embed- 
dings, we consider the spectrum of fluctuations along the Neumann-Dirichlet direc- 
tion a of the D5-brane as defined in (2.15) . Note that the term 'Neumann-Dirichlet' 
is valid only at large rh, where we can think of the AdS^ x 5*^ geometry as correspond- 
ing to the near-horizon limit of the gravitational background of a stack of coincident 
D3-branes since large m corresponds to a large 5*^ radius. To obtain the spectrum, we 
expand a = + (2TTa')6a in the DBI Lagrangian of the D5-brane and obtain the cor- 
responding quadratic Lagrangian. Here a is the polar angle on C AdS^ (2.15), with 
5*^ denoting the 3-sphere in the field theory directions. In general, fluctuations along a 
couple to fluctuations of the U{1) gauge field of the probe D5-brane (as discussed in 
detail in refs. [3, 4] for the "fiat" D3/D5 case). However, if one restricts the fluctuations 
to depend only on the time t and the holographic direction r, the modes decouple. It 
is the stability of this decoupled mode that we analyze. Another important point is 
that the corresponding gauge invariant operator is a dimension four (A = 4) operator. 
This is obtained from the asymptotic behaviour of the general solution of the equation 
of motion, using the standard field/operator dictionary, in analogy to the calculation 
we performed for the D7 brane case above. 

After implementing the ansatz 6a = e*'^*/i(r), solving numerically the corresponding 
equation of motion and imposing regularity of the solution, we obtain the plot of the 
spectrum versus bare mass presented in figure 8. 

At large rh, the spectrum should match the spectrum of the "flat" D3/D5 case. 
For the decoupled mode subject to our investigation, the analytic result obtained in 
refs. [3, 4] is given by: 

2777, / 

w = ^V(n + 3/2)(n + 5/2), (3.12) 

represented by the dashed black lines in figure 8. One can see the good agreement at 
large rh. Near the phase transition, the spectrum remains tachyon-free. At the phase 
transition it has a finite jump, as shown in the second plot in figure 8. Interestingly, 
at zero bare mass (m = 0) the spectrum is again described by equation (3.10) with 
A = 4. Once again the energy of the ground state is determined by its engineering 
dimension. A detailed analytic study of the spectrum at zero bare quark mass is 



-17- 




I 
I 
I 



8 










11=2 














7 
6 










n=l 






5 










11=0 








1.196 


1.198 


'1.200 


1.202 


1.204 



ill 



Figure 8: A plot of the spectrum Cj versus the bare quark mass rh for the spectrum of 
fluctuations corresponding to a. At large m the spectrum is described by equation (3.12), 
while at zero bare mass (m = 0) the spectrum is equidistant and is described by equation 
(3.10) with A = 4. The spectrum is tachyon free and has a finite jump at the phase transition. 

provided in sections 4 and 5. The lack of tachyon modes in the spectrum suggests that 
the D5-brane embedding is stable and does not slip off the equatorial two sphere of 
the three sphere where the field theory is defined (a = 0, S"^ C S^). 

4. Spectrum at zero bare mass — Gravity side 

Here we provide an analytic derivation of the spectrum of fluctuations at zero bare 
mass on the gravity side of the correspondence. To obtain the spectrum we recast 
the second order linear differential equation of motion of the fluctuation modes to one 
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dimensional Schrodinger equation. Our analysis is along the lines of the one employed 
in ref. [8] for the study of quasi-normal modes in the meson melting phase transition 
in flavoured supersymmetric Yang-Mills plasma. Let us briefly describe the method. 
We start with a second order differential equation of the form: 

h"{~p) + C,{~p)h\p) + {B{pfu' + Co(p))Mp) = . (4.1) 

Upon the substitution h{p) = a{p)f{p), with cr(p) defined via 

^^C,(p) + ^V (4.2) 



a(p) 2\ B(p) 

equation (4.1) becomes: 



^ ^ ^ "^-JiP) ] + {^'- Kfr(p)) /(p) = , (4.3) 



B{p) dp \B{p) dp 
where the effective potential K;fr(p) is given by 

K»(a = -^(c„(p) + C,(p)^ + ^) . (4.4) 

The final step is to define a new variable u{p) satisfying dz/dp = B{p). The resulting 
equation of motion is equivalent to a one dimensional Schrodinger equation for a particle 
with 2'mE = o)^ in the effective potential Vcs{z): 

f\z) + [uj' - V,^{z)]f\z) = . (4.5) 

This approach is particularly useful since it enables us to study properties of the spec- 
trum, such as the existence of bound states, even without solving the equation of 
motion. In our case the corresponding effective potential is relatively simple and we 
are able to obtain the spectrum in closed form. Let us present our results for the 
different modes studied in section 3. 



4.1 Fluctuations of a D7— brane probe 
4.1.1 Fluctuations along L 

We consider the ansatz SL = e*s*/i(p)3)' (5*^)3^' (S*^) in the equation of motion derived 
from the quadratic lagrangian (3.2). Here p = p/R, while 3^' and 3^' are spherical 
harmonics on C AdSs and C respectively. The corresponding differential 
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equation for h{p) is of the general form (4.1) with coefficients: 

^ l{i+2) _ l{l + 2) 8 + 16p^ 

5 + 16p^ + 48p^ 4 

-p+i6p^ ' ^^^) = rT^- ^'-^^ 

Applying the prescription outlined in equations (4.1)-(4.5) results in the effective po- 
tential 

3 + 4/(/ + 2) 3 + 4/(/'+2) , , , 

4 sm^ z 4 cos^ 2; 

where 

z(p) = 2arctan(2p). (4.8) 

The effective potential diverges at the interval boundaries z = 7r/2,7r, implying that 
the corresponding "wave function" f{z) should vanish there. It turns out that for this 
effective potential, we can solve equation (4.5) exactly. The solution regular at 2; = 7r/2 
is given by: 

f^{z) = C(cosz)'+52Fi[^(3 + [+/-w),^(3 + /~+/ + cj),/ + 2,cos2z](sinz)'+i . (4.9) 

Regularity at z = tt requires that one of the first two arguments of the hypergeometric 
function is a non-positive integer. Without loss of generality, we consider only positive 
Co. This implies that 

^(3 + / + /-w) = -n; n = 0,l,2,.... (4.10) 

Therefore we obtain the following expression for the spectrum: 

u = 3 + l + 2n + !] n = 0,l,2, .... (4.11) 

At / = / = 0, (4.11) agrees with the numerical analysis considered in section 3 above. 1 



4.1.2 Fluctuations along </>. 

In this subsection we study the spectrum of fiuctuations along (p in the limit of zero 
bare quark mass. Our starting point is equation (3.5). In this case there is a subtlety 
in taking the limit Z — in equation (3.5). Indeed if we substitute Z(p) = e^(p) into 
equation (3.5) and take the limit e — > 0, it is easy to verify that the coefficients of 
the equation of motion will remain ^(p) dependent. Indeed if we bring the equation of 



- 20 - 



motion to the general form (4.1), with dimensionless variables uj = uR and p = p/R), 
the coefficients are: 

The resulting effective potential Vcs{p) defined via equation (4.4) depends on the func- 
tion ^{p) and its derivatives. However C,{p) should also satisfy the linearized equation 
of motion for the classical D7-brane embedding. The next step is thus to solve for the 
highest derivative of ^ (p) from the linearized classical equation of motion and substitute 
into the expression for the effective potential V^fr(p). It is a straightforward exercise to 
verify that this completely removes the ^(p) dependence of the effective potential. 

Note that the function B[p) is the same as in equation (4.6) and hence z{p) is given 
by equation (4.8). In this coordinate z, the equation of motion is of the general form 
(4.5) with Vca{z) given by equation (4.7). Therefore we conclude that the spectrum of 
fluctuations along at zero bare mass is identical to the spectrum of fluctuations along 
L and is given by equation (4.11) which we duplicate below: 

u = 3 + l + 2n + l] n = 0,l,2, .... (4.13) 

One can see that at / = equation (4.13) agrees with the numerical results for the 
spectrum of fluctuations along from Figure 5 in Section 3. 



4.1.3 Fluctuations of the gauge field. 

In this subsection we are interested in studying the fluctuations of the U{1) gauge field 
of the probe D7-brane at zero bare mass. Our starting point is the action of the probe 
D7-brane: 

S=fJ deV\Gab + {27ra'yF,,\ + ^^^pr J P[C^,)] A F(2) A F(2) , (4.14) 
where -P[C4] is the pullback of the R-R four form whose "electric part" is defined by 

ciC(4) = ^VoKAdSs) , (4.15) 

where Vol(AdS5) is the volume form of AdSs. Note that in (4.14) we have written 
only the quadratic contribution (in a') to the Wess-Zumino term of the action. The 



- 21 - 



equation of motion is: 



^V.F- - Ip.^ (l + ^) (l - (l + ^^') V¥Ae'-d^A, = , (4.16) 

where a, 6 are general indices for the eight worldvolume coordinates of the probe D7- 
brane, while j, k are indices for the coordinates parametrizing the C 5*^ wrapped 
by the D7-brane. Note also that the Levi-Civita symbol e^^'^ is a tensor density (takes 
values 0, ±1) and ^3 is the determinant of the metric on the unit S'^ C AdSs. In 
this section we are interested in the spectrum of fluctuations at zero bare quark mass 
(L(p) = 0). Thus we write 

v^ = i?'^(p)v^v^; ^(^^) = ^^'(i + ^) (1-^) ' ^ = |- (4-1'^) 

Here is the determinant of the metric on the unit C . With these notations 
the equation of motion (4.16) can be written as 

5a(^(p) V^v^F"') - ^^g{P)^Ae'''^d,A, = . (4.18) 

Since we are in global coordinates, our probe brane is wrapping two three-spheres: one 
in the AdSs part of the geometry, C AdSs, and one in the part of the geometry, 
G S^. Moreover, we have to expand both in scalar and in vector spherical harmonics 
on the S^. Let us briefly review some basic properties of spherical harmonics: 



The scalar spherical harmonics y^Q^) on satisfy: 

ViV^y = -/(/ + 2)y , (4.19) 

and form a (|, |) irreducible representation of SO (A). On the other hand, there are 
three types of vector spherical harmonics on 5''^ [2]: The first type is constructed by 
simply applying a derivative on the scalar spherical harmonics d^y^ (longitudinal vector 
spherical harmonics) and transform in the (|, |) irreducible representation of 5*0(4). 
The other two types yl'^^fl^) (transverse vector spherical harmonics) transform in the 
(^, ^) irreducible representations and satisfy: 

ViV^y}^ - npl;^ = -(/ + ify-'^ , (4.20) 
^.'■'V,3;f = ±(/ + 1)^'^ , (4.21) 
V'yl'^ = . (4.22) 
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Here TZ^ = 251 Ricci scalar curvature of an unit S"^. Next we classify the fluctu- 

ations of the U{1) gauge field into three different types: 

Type I comprises of fluctuations Af along the 5*^ C which have expansion in 
transverse vector spherical harmonics 3^-'^. Since these modes satisfy VM^ = 0, they 
decouple from the remaining modes. Hence we have the ansatz: 

Af = ^f{p)e'^'y\~S'')y^^{S^)\ ^0 = 0; A„ = 0; A, = . (4.23) 

Here greek indices denote components along the 5*^ C AdSs. 

Type II are fluctuations of the gauge field along the C AdSs, which have 
an expansion in transverse vector spherical harmonics j^j-'^. Similarly to type I modes, 
these modes satisfy V"v4q, = and decouple from the rest. This implies the ansatz: 

A^ = ^f,{p)e''^'y]^CS'')y\S^); Ao = 0; A, = 0; = . (4.24) 

Finally type III modes are modes that do not fall into type I or type II classes. One 
may expect that in this case all modes are coupled and should be studied simultaneously. 
However only vector spherical harmonics transforming as scalar spherical harmonics 
can be used, which are given by dcy\S'^) and by diy\S^). Therefore both A^ and 
Ai components of the gauge field are pure gauge and one of these sets may always be 
gauged away. Thus without loss of generality, the type HI modes can be described by 
the ansatz: 

Ao = ^m{p)e''''y\~S^)y\S^); A^ = <^mip)y\s')y{S'); (4.25) 

A = ^ni{p)y\s')diy{s'y, a^ = o. 

Next we proceed with the analysis of the equations of motion. Let us begin with 
type I modes. After substituting the ansatz (4.23) into the equations of motion (4.18) 
we obtain: 



{Gipyp^^f) + 



cl;2 l{l + 2) (/ + 1)(/ + 1±4) 



ip' + \? ip' - ly p' 



<l>f = 0. 

(4.26) 



Equation (4.26) is of the general form (4.1) with coefficients given by: 



3 + 16p^ + 80p^ 4 l{l + 2) (/ + 1)(/ + 1±4) 

^'^P^ = pW-1) ' = VTI ^'^P^ = ~W^' P ■ 

(4.27) 
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Using the procedure outlined in equations (4.1)-(4.5), we obtain an equation of the 
form (4.5) with effective potential 

^ , (2/ + l±4)(2/ + 3±4) 3 + 4/7/ + 2) , , , ^ 

where z{p) is given by equation (4.8). The general solution regular at 2 = 7r/2 is given 
by: 



/±(^) = C±(cosz) 



, „ 7 , 1 



(4.29) 



X 2i^i[-(3 ± 2 + / + / - cj), -(3 ± 2 + / + / + w), / + 2, cos^ z](sinz)'^2 



Regularity at 2 = vr requires that one of the first two arguments of the hypergeometric 
function must be a non-positive integer. Without loss of generality we can consider 
only positive a) and find 

-(3±2 + /~+/-c:;) = -n; n = 0,1,2,.... (4.30) 

Therefore we obtain the following expression for the spectrum: 

= 3±2 + / + 2n + [; n = 0,1,2,.... (4.31) 

At large p the equations of motion (4.26) are the same as the equations of motion of 
type I modes for the "flat case" studied in ref. [2]. Therefore the conformal dimension 
of the corresponding gauge invariant operators are the same, namely A± = 3 ± 2 + /. 
Finally for the spectrum of type I fluctuations we can write: 

cu^ = A± + 2n + / n = 0,1,2,... . (4.32) 

We see that similarly to the spectrum of fluctuations of the transverse scalars, the 
ground state is again determined by the conformal dimension of the corresponding 
gauge invariant operators. In the next section we will confirm equation (4.32) from 
field theory considerations. 



We now focus on the spectrum of type II modes. Substitution of the ansatz (4.24) 



into the equations of motion (4.18) results in the following equation: 



16p 



(16p4 - 1 



-dr, 



(16p^-i; 
16p 



([+1)2 l{l + 2) 



p2 



= 



(4.33) 



which is of the type (4.1) with coefficients: 

1 + 48p^ 



Clip) 



p — 16p^ 



B{p) 



4p2 + 1 



Co(p) 



([+1)2 l{l + 2) 



(4p 



n2 



p2 



(4.34) 



Following the procedure outlined in (4.1)-(4.5) one obtains an equation of the form 
(4.5) with effective potential 



Kff(^) 



3 + 4/'(/~+2) 3 + 4/(^ + 2) 



4 cos^ z 



4 sin^ z 



(4.35) 



where z{p) is given by equation (4.8). The general solution regular at 2 = 7r/2 is given 
by: 

/ [z] = C(cos z) i+'aFi (3 + /" + / - w) , ^ (3 + [ + / + w) , / + 2, cos^ z] (sin 2) . (4.36) 

Regularity at z = vr requires that one of the first two arguments of the hypergeometric 
function be a non-positive integer. Therefore the spectrum of type II modes is given 
by: 

u = 3 + l + 2n + l] n = 0, 1, 2, . . . . (4.37) 

Given that the conformal dimension of the dual gauge invariant operator is A = / + 3 
[2] one obtains: 



w = A + 2n + /; 



n = 0,1,2,... . 



(4.38) 



Once again the ground state is determined by the conformal dimension of the dual 
gauge invariant operator. 



Finally we study the spectrum of type III modes. Substitution of the ansatz (4.25) 
into (4.18) results in the following equations coming from components of the gauge field 



along the t, p, C and C AdSs directions: 



'-^2^1,2 '^in + ^^^p-^m) - f~2 iA2 '^"i ^("^in - iii) 



4 



^^a.(^*.„) - ^4„. - ^*„, ^ . (4,39) 

P^G{p) - i)2 (p2 + 1)2 p2 



Note that the third and fourth equations in (4.39) are vahd only for / > and / > 0, 
respectively. Moreover, note that we have four equations for three unknown functions 
'^'iiij '^'iii and $111. However one can show that only three of the equations are indepen- 
dent. In particular, it is easy to verify with simple algebraic manipulations that the last 
three equations in (4.39) imply the first one, which may thus be skipped. Moreover, 
generically u ^ and one can use the fourth equation to eliminate <I>ni from the second 
and third equations in (4.39). The resulting system of coupled equations for $ni and 
$111 is given by: 

Klip) + CnipWmip) + [BiipyCo^ + Coi(p)]$iii + Mu{p)^m{p) = , 
Klip) + Ci2{p)Kiip) + [52(p)'cu2 + Co2(p)]$iii + M2i(p)$ni(p) = , (4.40) 

where the coefficient functions are given by: 

1 - 16p^ + U2p' 4 2/(/ + 2)(l-4p^) 

Cuip) = ^ _ 16^-5 ^ ^^(^) = rTv' = — WTW — ' 

r (~.^- ^ + ^^P' + (1 + 12p^ + 72p") Kl + 2) Kl + 2) 

(p - 16p5)2 (p2 _ 1)2 p2 

3 + 16p^ + 80p^ 4 (2 + 8p2) 

<^i2(p) = T — ; B2{p) = — — TTj; M2i(p) 



C^02(P) 



p-16p5 ' l + 4p2' p-Ap^ 

1(1+2) l{l + 2) 



ip' - P' 



Note that the functions Bi and B2 in equations (4.40) are the same. This suggests that 
we may employ similar procedure to the one outlined in equations (4.1)-(4.5). To this 
end we define: 

*m(p) = /i(p)(Ti(p); $m(p) = f2{p)(T2{p) , (4.42) 
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where ai and ct2 are defined via: 



The system (4.40) may be brought to the form 



-Ucu + §]; 1 = 1,2-. (4.43) 



f^{z) - V}i\z)]f,{z) + rm,{z)Uz) = , (4.44) 
f^{z) +[Co' - vS\z)]Uz) + m,,{z)f,{z) = 0; 



where 



rn,.i~~)--m^; (4.45) 

4sm"'2; 4cos2 2; sm^z 



4 sm 2 4 cos^ z sm 2; 

and z{p) is given by equation (4.8). Note that the difference of the effective potentials 
AV^ff = vjg — V^g^ and the "coupling functions" mu, 77121 are all proportional to 
1 / sin^ z and differ only by constant multiplicative factors. This fortunate property 
enable us to decouple easily the system of equations (4.45). Indeed let us form the linear 
combination ?7 = /i + a/2 for some constant parameter a. Combining the equations in 
(4.45) results in the following equation of motion for ri{z): 

V"iz) + [uj' - V^^^Uz) + am,,{h + !!^ii^^^/,) = . (4.46) 

am2i 

Note that due to the property mentioned above, the fraction (mi2 — aAVi3fr)/(am2i) is 
a constant and therefore we can demand that it is equal to a. We may now write 

ri"{z) + [Q^ - V,s{z)]ti{z) = 0; V,^{z) = V^^\z) - am2i{z) . (4.47) 

The parameter a is determined by: 
(mi2 - aAVeflf) 



(am2i) 



^ - 2a - /(/ + 2) = 0; ^ a+ = / + 2; a_ = -/; . (4.48) 



Therefore we end up with two decoupled modes 'q±{z) = fi{z) + a±f2(z) described by 
the effective potentials: 

. + i±«i±^ . (4.49, 

''^ 4sm^z 4cos2 2; 
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The corresponding solutions of the equations of motion regular at 7r/2 are given by: 

X2Fi[^(3 ±1 + 1 + 1- ^(3 ± 1 + 1 + I + c^±), / + 2, cos^ z]{sm z)'^^^+^ . 

(4.50) 

Regularity of the solutions at ^ = vr implies the following eigenfrequencies: 

d)± = 3 ± 1 + / + 2n± + /; n± = 0, 1, 2, . . . ; . (4.51) 

Note that the spectrum of ?7+ is contained in the spectrum of the ?7_ modes. Therefore 
modulo the number of independent modes the spectrum of type III modes is described 
by: 

u = 2 + 1 + 2n + l, n = 0,l,2, . (4.52) 

On the other hand, by studying the asymptotic behaviour of the solutions at large p, 
we can verify that the conformal dimension of the dual gauge invariant operators is 
A = 3 + /. Therefore we can write our expression for the spectrum as: 

u = A-l + 2n + l. (4.53) 

Note that the entire spectrum is shifted by —1 relative to the spectrum of type II 
modes. The spectrum of modes with / = and / = requires separate analysis. Here 
we will simply report the result, which is that the expression for the spectrum is still 
given by equation (4.52) at vanishing / or /. We will comment more on the spectrum 
of type III modes in the Section 5 below, where we reproduce our results from field 
theory considerations. 

4.2 Fluctuations of a D5— brane probe 

In this subsection we obtain the spectrum of fluctuations of a probe D5-brane at 
vanishing bare mass m = in closed form. Let us first consider fluctuations of l{r). 

4.2.1 Fluctuations of l{r) 

We consider the ansatz 61 = e^Ti^h{f)y\S'^)y\S'^) for the fluctuation of the D5-brane 
and substitute this into the corresponding equation of motion (here r = r/R) . The 
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result is a differential equation of the general form (4.1) with coefficients: 



Co(p) 



H J f-16f5 ' ^ ^ l + 4f2 ^ ^ 

The calculation proceeds in complete analogy to the D3/D7 case. One applies the 
procedure outlined in equations (4.1)-(4.5) and obtains the effective potential: 

V^siz) = + z e [7r/2, vr) , (4.55) 

sm z cos^ z 

where 

z{r) = 2 arctan (2f). (4.56) 
The solution regular at ^ = 7r/2 is given by: 

f^{z) = C(cos zy+\F^[^{2 + !+l-Q),^{2 + l + l + Q), 1 + ^, cos^ z](sin zY^' (4.57) 

Regularity at z = vr requires that one of the first two arguments of the hypergeometric 
function must be a non-positive integer. Without loss of generality we can consider 
only positive u, which implies 

-{2 + 1 + 1 -u) = -n; n = 0,1,2,... . (4.58) 

Therefore we obtain the following expression for the spectrum: 

u = 2 + l + 2n + !] n = 0, 1, 2, . . . . (4.59) 

Equation (4.59) is equivalent to equation (3.10) with A = 2, / = 0, / = 0. In section 5 
we will obtain this equation from field theory considerations. Now let us focus on the 
fluctuation of the Neumann-Dirichlet coordinate. 



4.2.2 Fluctuation of the Neumann-Dirichlet coordinate. 

In this subsection we analyze the spectrum of fluctuations along the polar coordinate 
a defined in (2.15) of the three-sphere inside the global AdSs space-time. Note that 
in the limit of infinite radius of the three-sphere, we obtain AdSs space-time in a 
Poincare patch. In this limit we have an underlying description in terms of a stack 
of Nc coincident D3-branes and the notion of Neumann-Dirichlet coordinate is well- 
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defined. 



To obtain the spectrum of fluctuations we follow tlie same prescription tliat we did 
so far. We substitute 5a = e^^h{f) into tlie equation of motion for Sa (here f = r/R) 
and obtain a differential equation of the general form (4.1) with coefficients: 



Coif) 



32 



4 + 24f^ + 96f^ 
Ci{r) = —; B[r) 



;i-4f2)2' ^^^'^ f(16f4-l) ' ^' ' l + 4f2 

The corresponding effective potential Ve^if) defined in equation (4.4) is given by 

3(1 +4f=2^2 



(4.60) 



Kff(r) 



(4.61) 



After changing variables to u{r) = 2 arctan (2r) the equation of motion takes the general 
form (4.5) with effective potential Kff(w) given by 



The general solution is given by: 



6 



• 2 ' 

sm u 



u e {n/2,n) . 



f{u) = DiPI'\^^{cosu) + L>2<5^„\/2(cos 



u 



a/sI 



smw 



(4.62) 



(4.63) 



In the original variables {h{p)), we can verify that regularity of the solution requires 
that /(7r/2) = /(tt) = 0. This requires: 



cos 




TT 

2 


sin 








COS 




TT 

2 


sin 






D2_ 







(4.64) 



Non-trivial solutions exist only if the determinant in (4.64) vanishes. This is equivalent 
to sin(a)|) = and hence Co = 2ni, where ni = 0,1,2, ... .In fact, even with this 
restriction we can verify that h{f) diverges at f = 1/2. In particular, it behaves as 
oc l/(f — 1/2). However this still corresponds to bounded fluctuation of the probe D5- 



brane since the invariant quantity Gaai^a) 



oc r 



l/2)2/i(r)2 < (yo remains bounded. 



Next we impose the condition that at infinity (f — > oo) h{r) falls as oc 1/f^. This 
is related to the fact that the corresponding gauge invariant operator has conformal 
dimension four (A = 4). One can verify that all solutions with ni > 2 have that 
asymptotic behaviour. Therefore the final expression for the spectrum is: 



n 



0,1,2, 



(4.65) 
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which is exactly equation (3.10) with (A = 4). 



5. Spectrum at zero bare mass — Field theory side 

We now turn to the analysis of the spectrum at zero bare quark mass on the quantum 
field theory side of the correspondence. It is instructive to begin with a brief review of 
the properties of a conformally coupled free scalar on x S*". 

5.1 A Conformally coupled scalar on M} x S"" 

We begin by reviewing the derivation of the spectrum of a conformally coupled free 
scalar field on an Einstein universe M} x S"^. We start with the Lagrangian 



l9l!:)da<i>d,<i> + Un^ 



(5.1) 



Here g{n)af, is the metric on M} x S*" given by: 

ds"" = -df + Rldnl , (5.2) 
TZ^^^ is the Ricci scalar curvature of the metric in equation (5.2) given by: 

n(-) = n{n-l)-^ (5.3) 

and Rn is the radius of 5*". For ^ = {n — l)/4(n), the scalar is conformally coupled. 
Now it is straightforward to write down the equation of motion: 

-^00 + - ^-^152^ = • (5-4) 
Next we Fourier transform : 

(j){t,nn)= f dujJ2Ci{oo)e''''y{^ny, / = 0,1,...; (5.5) 

I 

where 3^'(f2„) are the scalar spherical harmonics on the unit 5"" satisfying: 

An,yinn) = -lil + n-l)y{Qny, /~=0,1,...; . (5.6) 
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The equation of motion (5.4) becomes: 



cj^ - f(n - l)/2 + f ^ ^ 



-"'n J 



= , (5.7) 



which imphes: 
where Ej is given by: 



Q(w) oc 6{uj' - Ef) , (5. 



Ei= ((n-l)/2 + /) /'=0,1,...; . (5.9) 



Now one can write equation (5.5) as: 

- \ / J--^T 
I V I 

and proceed with quantization of in the standard way. 
5.2 A collection of fields on 5"^ 

In this section we calculate the spectrum of a collection of free scalar fields propagating 
on S^. In particular we consider the field content of the gauge invariant operators holo- 
graphically dual to the normal modes of the D7-brane that we computed in section 4. 
The flavoured dual field theory in question is a A/" = 4 SYM theory in 1 + 3-dimensions 
coupled to an A/" = 2 hypermultiplet. The coupling of theory is conveniently studied 
if we decompose the A/" = 4 supermultiplet of the adjoint degrees of freedom into one 
vector multiplet and two M = 2 hypermultiplets. The bosonic content of the A/" = 4 
multiplet is the (1 + 3)-dimensional vector and the six adjoint scalars X^, . . . X'^ . 
The bosonic sector of the vector multiplet contains the vector and the adjoint scalars 
Xy = (X^, X^). The rest four adjoint scalars X^ {i = 4, 5, 6, 7) form the bosonic con- 
tent of the two adjoint hypermultiplets. The fundamental fields form a hypermultiplet 
comprised of two spinor fields ipi and two complex scalar fields g"^ {i,m = 1,2). 

The operators corresponding to mesonic excitations of flavoured M = 2 SYM theory 
on flat Minkowski spacetime have been identified in refs. [2], [22]. In particular the 
meson states have been classified in multiplets forming irreducible representations of 
the global 5*0(4) S'f/$(2) x SUr{2) symmetry. Furthermore it has been verified that 
states from the same supermultiplet (sharing the same charge under SUr{2) ) have the 
same mass spectrum. The latter has beeg determined via holographic techniques. 

The flavoured SYM theory on M^x that we consider has the same global SO {4) ^ 
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SU^{2) X SUr{2) symmetry. This is why we may use the same classification of the 
mesonic excitations. However, the supersymmetry algebra on X 5* mixes states at 
different floors of the Kaluza-Klein tower on 5*^. Therefore we cannot expect to have 
the same degeneracy of the spectrum of meson excitations as for the theory on flat 
Minkowski spacetime. A summary of our results is provided in table 1. To begin with, 





fluctuation 


d.o.f. 




op. 


A 


spectrum 


mesons 


1 scalar 


1 


V 2 ' 2 ' 




1+i 


A + 2n + i 


(bosons) 


2 scalars 


2 


V2' 2^ 


Mf 


3 + i 


A + 2n + i 




1 scalar 


1 


V2' 2^ 




3 + i 


A + l + 2n + £ 




1 scalar 


1 


{- -) 
V2' 21 


jO£ 


3 + £ 


A-l + 2n + i 




1 tr. vector on 


2 


{- -) 
V2' 21 




3 + i 


A + 2n + i 




1 scalar 


1 


(£+1 

V 2 ' 2 / 




5 + i 


A + 2n + i 



Table 1: Bosonic sector of mesonic excitations of conformally coupled flavoured J\f = 2 
SYM on X 5^ . Here {j^,jR) labels irreducible representations of the global 50(4) 
SU^{2) X SUji{2) symmetry, while i labels the multipole momentum on . Note also that 
the transverse vector satisfies ^ aJ'^^ = (Coulomb gauge), where Vq, is a covariant 
derivative on . The hyphen in the fifth column corresponds to a descendent of the operator 
C^^. Note also that the spectrum of the pseudo-scalar J^^ and that of the descendent of C^^ 
has been obtained only from supergravity in this paper, while all others are obtained on both 
sides of the correspondence. 



let us consider the operator dual to fluctuations of the transverse scalar modes of the 
probe brane (fluctuations along L and 0). The corresponding operator is given by [22]: 

= A(Tf,x'^j + qmX^x'qm (2, m = i, 2) , (5.11) 

which has conformal dimension A = / + 3. Here denotes the symmetric traceless 
operator insertion X^^^ ...X"^^^ of / adjoint scalars (z = 4,5,6,7), Xy denotes the 
vector (X*,X^) and = ((7^,cr^) is a doublet of Pauli matrices. Let us focus on the 
scalar part of the operator (5.11). Without loss of generality we can consider only the 
{m = 1, A = 8) additive term. In view of the result obtained in the previous subsection 
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we expand: 

X^Pit) = \m[e^^X^., + e-"^^X^\): C^ = 4,5,6,7,8; (5.12) 
[1=0 h=o V2 + 2/i 

l,=0 l2=0 V2 + 2/2 

Here we have written explicitly the i,j = 1 . . . Nf. color indices. Note that we have 
turned on only the "zero" mode of the adjoint scalar which is due to the conformal 
couphng. Next we define a vacuum state |0) satisfying: 

= blrjO) = ^^\0) = . (5.13) 

The state representing the collection of fields corresponding to the {m = 1, A = 7) 
component of the operator (5.11) is given by 

[1/2=0 ^1=0 h=o ^J {I2 + l)([i + 1) 

(5.14) 

Note that the state (5.14) is a superposition of states with definite energy Ei j = 
(3 + / + J)/R3, where J = h + I2. Apparently there is a large degeneracy of the 
spectrum corresponding to the different choices of li and I2 which sum up to the same 
number J. Note also that a state with a definite energy has a definite eigenvalue under 
the antipodal map (which for spheres in even dimensions coincides with parity) on the 
three-sphere given by (—1)"^. Our next step is to expand the product of the spherical 
harmonics in (5.14) in Laplace series: 

00 

y./.(^3)yV.(^3) = ^ ^ dl^^^^,^^/\S') . (5.15) 

[=0 ^=0 

The coefficients C'/, r , are non-zero only for IL — lol < I < h + h (addition of angular 

momentum) and (—1)' = (—1)'^"^'^ (conservation of the antipodal map eigenvalue)^. 
Therefore the Laplace series in (5.15) terminates at J = /i + 12- For all values of / that 
^We refer the reader to the appendix of ref. [23] for a detailed review of the properties of SU (2) x 



- 34 - 



appear in the expansion, we can write J = 2n + l for some non-negative integer n. This 
imphes that a state with a definite energy E j^i can be expanded as: 

\Ej,) = l^e^^-* E y'(S')Cl,^\0) , (5.16) 

2n+l=J 1=0 

where C^- is defined by: 

nil 

Ol + l)" ih + lf Qll 

E E E . 4,.<X^^'-l,.. (5.17) 

[l+[2=2n+[ ^1=0 h=0 y (/2 + l)(/l + 1) 

Now the general state (5.14) can be written as: 

qiXW^ilO) = f E E e^^^^'^'"^'^^y"^(5^)C^LlO) • (5-18) 

n,l=0 ^=0 

Therefore if we define u = uR^, our final expression for the dispersion relation is given 
by 

u = 3 + l + 2n + l; n = 0,l,2, .... (5.19) 

This coincides precisely with equation (4.11) which is obtained for the dual mode on 
the gravity side of the correspondence. The fact that our free field theory analysis 
reproduces the spectrum at strong coupling is consistent with the non-renormalization 
theorem for chiral primaries. Note that this matching provides a non-trivial check of 
the holographic technique employed in Section 4. 

Vector modes. Next we apply the same approach to study the spectrum of gauge 
invariant operators corresponding to fiuctuations of the U{1) worldvolume gauge field 
of the D7-brane. There is only one chiral primary operator corresponding to the type 
I A~ mode analyzed in section 4. It has conformal dimension A_ = 1 + 1 and is given 
by^: 

= r^Lx'-'g" , (5.20) 

where cr^„ (/ = 1,2,3) is a triplet of Pauli matrices. To obtain the spectrum of the 
corresponding collection of scalar fields one may apply a procedure analogous to the 

SU{2) Clebsch-Gordan coefficients. Note that half integer quantum numbers rather than integer ones 
are used in ref. [23]. 

•^Note that we are using the convention for the quantum number I employed in ref. [2], namely the 
ground state has ^ = 1. 
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one outlined in equations (5. 12)- (5. 18). The ground state is again given by the sum of 
the zero point energies of the constituent fields, namely (l-|-l-|-(/ — l))/i?3 = (1 + Z)/-R3. 
Furthermore by expanding in Laplace series and employing properties of the Clebsch- 
Gordan coefficients on can argue that the spectrum is given by the negative sign in 
equation (4.32): 

uj- = A_ + 2n + l n = 0,l,2, .... (5.21) 

The operator corresponding to the Af mode is descendent of the operator C^' and can 
be analyzed in a similar way. 

Next we consider the spectrum of type II modes described in section 4. Note 
that our type II modes are subset of the type II modes defined in refs. [2] , [22] . The 
latter correspond to the components of the U (Nf) fiavour current operator [22] in 1 + 3 
dimensions: 

J^' = ^t^abXVi + ^QmX'D^qm ' tD^QmX'qm (/i = 0, 1, 2, 3) (5.22) 

and satisfy d^J^^ = (Lorentz gauge) which leaves only three independent degrees of 
freedom. In our case the field theory is defined on x and the Lorentz symmetry is 
broken. Our type II modes (look at section 4) comprise of the spacial part of the current 
operator (5.22) defined on and satisfy VaJ^°'^ = (Coulomb gauge). Therefore we 
have only two independent degrees of freedom which have Laplace expansion in terms 
of transverse vector spherical harmonics 3^^"^. It turns out that one can successfully 
apply the procedure outlined in (5. 12)- (5. 18) with minor modifications. Consider the 
scalar part of the operator (5.22) one can see that the energy levels of the corresponding 
collection of fields are given by (1 + 1 + / + /i + li)/Rs. The relevant Laplace expansion 
is given by: 

Note that the coefficients G^^ ^ ^ are non vanishing only for li + I2 = 2n + 1 + I [23] 
for non-negative integer values of n. Therefore the final expression for the spectrum of 
oj is given by equation (5.19) again. 

Finally we consider the spectrum of type III modes. The corresponding operators 
are the time component of the U{Nf) fiavour current (5.22) and the pseudo-scalar 
operator [22]: 

J''~' = ^hffex'-V- + ■ . • ; > 1) • (5.24) 

There are is also a contribution from the space-like components of the fiavour current 
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(5.22) which have expansion in longtudinal vector spherical harmonics Vq3^'(iS'^). The 
corresponding spectrum (modulo degeneracy of the ground state) was calculated in 
section 4 and is given by equation (4.52). To obtain the spectrum from field theory 
considerations, we consider the relevant contributions from the scalar fields to the 
flavour current (5.22). The energy levels of the corresponding collection of fields is given 
by {1 + 1 + 1 + li + li) / R-^. The relevant Laplace expansion is in terms of scalar spherical 
harmonics (for the time-like component of the flavour current) and longitudinal vector 
spherical harmonics Vq,3^' (for the space-like components of the flavour current). Using 
properties of the r , Clebsch-Gordan coefficients, we can show that in both cases 

hhhh 

the l-th terms in the expansion vanish unless li + I2 = 2n + I. Therefore the spectrum 
is given by 

w = 2 + / + 2n + / , (5.25) 

which coincides with the gravity result (4.52). Note that the type III vector modes 
analyzed in section 4 have two independent modes ri± with spectrum u}± given by 
equation (4.51): 

u± = 3±l + l + 2n + l, (5.26) 

and we verified that the contributions from the scalar fields to the flavour current (5.22) 
reproduce the spectrum of the r)- modes. It is natural to assume that the pseudo-scalar 
operator (5.24) corresponds to the 77+ modes and has spectrum given by the positive 
sign in equation (5.26). 

This completes our study of the meson spectrum of the bosonic sector of confor- 
mally coupled flavoured M = 2 SYM theory on x S''^. 

5.3 A collection of fields on S"^ 

In this subsection we consider the gauge field theory holographically dual to a D5- 
brane probe in a global AdSs x space-time. The theory in question is an AA = 4 
SYM theory on M} x coupled to a defect field theory defined on an equatorial 
S"^ C S^. At vanishing bare quark mass, the theory is superconformal and preserves 
A/" = 4 supersymmetry in 2+1 dimensions. The coupling of theory is the same as 
the one for the theory on flat Minkowski space studied for a first time in refs. [9, 10] 
and reviewed in ref. [24]. The analysis requires decomposing the (3 + l)-dimensional 
A" = 4 supermultiplet into two (2 + 1) -dimensional A/" = 4 supermultiplets, a vector 
multiplet and a hypermultiplet. The bosonic content of the (3 + l)-dimensional A/" = 4 
multiplet is the (3 + l)-dimensional vector and the six adjoint scalars X^, . . . X'^ . 
The bosonic content of the (2 + l)-dimensional vector multiplet comprises of the 2 + 1 
dimensional vector Ak k = 0,1,2 and the three scalars Xy = (X'^,A^,A^). The 
bosonic content of the (2+l)-dimensional hypermultiplet is the scalar A3 and the scalars 
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Xh = {X'^, X^, X^). The flavour fields form a (2 + l)-dimensional hypermultiplet with 
two fermions ipi and two complex scalars (i, m = 1, 2). 

A full analysis of the spectrum of mesonic excitations is beyond the scope of this 
paper. In this subsection we calculate the spectrum of fluctuations of the operator 
corresponding to the fundamental condensate of the defect field theory given by the 
operator [24]: 

= ^i4xHj + 2g™X^^x'g"; {A = 1, 2, 3; i,m = 1, 2) . (5.27) 

Here denotes the symmetric traceless operator insertion X^^"^ . . . X*'^ of / adjoint 
scalars {i = 4,5,6), Xy denotes the vector {X'^,X^,X^) and = {(r^,a'^,a^) is a 
triplet of Pauli matrices. On gravity side this operator is dual to the fluctuations 
modes of the transverse scalars of the probe D5-brane. For a particular component 
(say A = 7), the operator S^'' corresponds to fluctuations of the profile function /(r) 
studied in section 4. Therefore the corresponding spectrum is given by equation (4.59). 
Our goal is to obtain the same result from field theory considerations. Note that the 
operator S"^^ is the analog of the operator Aif- considered in the previous subsection. 
This implies that the spectrum can be obtained by adapting the procedure outlined in 
equations (5.12)-(5.18) to the defect field theory defined on the equatorial 5^. To begin 
with let us focus on the scalar sector of the A = 7 component of the operator E"^^ and 
consider only the m = 1 additive term. The analog of equation (5.12) is given by: 




(5.28 

(i=Omi=— <i ^ ^ 
,-_n™.=_/„ Vl + 2/2 



00 li 



Note that in equation (5.28) the adjoint fields X'^ are truncated to the zero modes 
on S^, while the fundamental scalars are expanded in spherical harmonics on S"^ and 
have zero point energies l/(2i?3). Note also that R2 = R3 since the fundamental fields 
propagate on an equatorial two-sphere S"^ inside the three-sphere S^. Next we define a 
vacuum state |0) satisfying: 

aU|0) = fcUlO)=^o%|0) = 0. (5.29) 
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The state representing the collection of fields corresponding to the {m = 1, A = 7) 
component of the operator (5.27) is given by: 

q,XW<li\0) = 

h,!2=0mi=-li m2=-l2 Y (2^2 + l)(2/l + 1) 

(5.30) 

Note that the state (5.30) is a superposition of states with definite energy Ei^j = 
(2 + / + J)/R3, where J = li + I2. Note also that a state with a definite energy has a 
definite parity given by (—1)'^. Next we expand: 

l=\l\—l2\ m=—l 

The coefficients Cr™ r are non-zero only for l/i— /2I < ^ < h + h (addition of angular 

Iimi,l2m2 — o 

momentum) and (—1)' = (—l)'^"*"'^ (conservation of parity). Therefore for all values of 
/ that appear in the expansion, we can write J = 2n + / for some non-negative integer 
n. This implies that a state with a definite energy Ej^i can be expanded as: 

2n+l=J m=—l 

where C\ is defined by: 

nlm 

CL- E E E ^ a^^X^x^H.. ■ (5-33) 

i-^+l2=2n+l mi=-h m.2=-l2 y (2^2 + l)(2/l + 1) 

Now the general state (5.30) can be written as: 

oo I 

q.XW^M = E E R.e^'^'^'''^'^^^y\~S')Clj^) . (5.34) 

n,l=0 m=—l 

Therefore if we define a) = uR^ , our final expression for the dispersion relation is given 
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precisely by equation (4.59): 

u = 2 + l + 2n + l] n = 0,l,2, .... (5.35) 

This completes our study of the meson spectrum of the defect field theory on S"^. The 
free field theory analysis of the spectrum agrees with the results obtained at strong 
coupling via holographic techniques. Similarly to the case of the flavoured gauge field 
theory on considered in the previous subsection, this is in agreement with a non- 
renormalization theorem. However the perfect matching of equations (4.59) and (5.35) 
provides a satisfying verification of the validity of the holographic approach towards 
strongly coupled systems at least in the presence of some supersymmetry. 



6. Conclusion 

We have considered meson spectra on R x S*^ in a holographic approach of embedding 
D7 or D5 brane probes into global AdS space. First, in a numerical approach, we de- 
termined the bare quark mass dependence of the spectra through the phase transition 
which occurs between branes which do or do not reach the 5*^ sphere in the field theory 
directions. Second, we studied the spectra at zero bare quark mass analytically and 
found exact agreement between the gravity and field theory calculations. For the field 
theory calculation, we expanded the fields in the operators considered on and com- 
bined them using the relevant Clebsch-Gordan coefficients together with a generalized 
parity argument. The fact that the free field calculation agrees with the gravity calcu- 
lation is consistent with the non-renormalization theorems of supersymmetric theories. 
Moreover, our calculation confirms the physical interpretation that since the zero-point 
energy on exceeds the binding energy of the mesons, the mesons cannot pair-produce 
any more. 

As far as we know, this is the first example of a calculation within gauge/gravity 
duality with added flavour using the Born-Infeld action where there is an exact match 
between a non-trivial field theory calculation and the gravity result. In the plane wave 
limit, field theories similar to ours have been studied previously in [25, 26, 27]. We 
hope that the results of this paper will inspire further field-theory calculations in top- 
down gauge/gravity models where the dual field theory is known explicitly, in view 
of further comparisons between weak and strong coupling results for observables of 
physical relevance. 
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